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Abstract 

A quantum system is described, whose wave function has a complexity which increases 
exponentially with time. Namely, for any fixed orthonormal basis, the number of compo- 
nents required for an accurate representation of the wave function increases exponentially. 
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This Letter describes a quantum system whose representation becomes increasingly 
complex with the passage of time. This behavior, which is well known as a generic 
property of classical Hamiltonian systems, is commonly called "chaos." For example, an 
initially smooth Liouville density becomes more and more convoluted as time passes, and 
it may form intricate shapes with exceedingly thin and long protuberances. However, 
quantum wave functions usually have a milder behavior [1]. In particular, quantum 
systems with discrete spectra can be represented, with arbitrary accuracy, by a finite 
number of eigenfunctions which are periodic in time. Their wave function thus is almost 
periodic. Its computational complexity [2, 3] does not increase as time passes. (Some 
authors use other criteria for complexity, such as the visual appearance of wave functions, 
but these seem rather subjective.) 

In this Letter, I describe a quantum system whose wave function becomes increasingly 
complex with time, just as a Liouville density would do in classical mechanics. The long 
range evolution of the wave function is effectively unpredictable with finite computing 
resources. This is a genuine example of quantum chaos. I have no pretense that the 
abstract model discussed here represents, even approximately, a real physical object. 
However, this is a formal proof that, contrary to some claims, quantum mechanics is 
compatible with the existence of chaos. 

Such a chaotic quantum system can be constructed for any classical area-preserving 
map. As a concrete example, I shall choose a map consisting of alternating twists and 
turns of a unit sphere. That map has a quantum version, called "kicked top" [4], with 
a discrete, finite-dimensional quasi-energy spectrum. Therefore, unlike the classical map, 
the quantum kicked top is not truly chaotic. However, a classical dynamical system may 
have several, different quantum versions. Another quantization of the same twist and 
turn (T&T) map is presented here, which is just as chaotic as the original classical map. 

The latter is defined as follows: consider the unit sphere, x 2 + y 2 + z 2 = 1. Each step of 
the T&T map consists of a twist by an angle a around the z-axis (namely, every xy plane 
turns by an angle az), followed by a 90° rigid rotation around the y-axis. The result is 



z' = —x cos(az) + y sm(az). 

This map is obviously area preserving. For low values of a, most classical orbits are 
regular (that is, they are quasi-periodic). As a increases, so does the fraction of chaotic 
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orbits, until for a = 3 most of the sphere is visited by a single chaotic orbit (all numerical 
calculations below refer to the case a = 3). 

The "kicked top" [4] is a mechanical system (classical or quantal) which mimics the 
above geometrical map. In the classical case, the top has a 3-dimensional generalized phase 
space [5], with canonical variables J x , J y , and J z , satisfying Poisson brackets { J x , J y } = 
J z , and cyclic permutations. The mapping 

V = T 

J' y = J x sin(aJ 2 /J) + J y cos(aJ z /J), (2) 
J' z = —J x cos(aJ z /J) + J y sin(aJ 2 /J), 

is a canonical transformation [6] which leaves J 2 invariant. The classical values of J x /J, 
Jy/J, and J z /J, lie on a unit sphere, and therefore they transform just as x, y and z in 
Eq. (1). 

A natural way of quantizing these equations is to replace the classical variables Jk 
by the corresponding quantum operators. Equation (2) then becomes a quantum map, 
generated by the unitary transformation 

U = exp{-i7Tj y /2h) exp{-iaJ 2 /2Jh). (3) 

For any eigenstate of J, this U is a matrix of order (2j + 1), with a discrete spectrum. 
The evolution is multiply periodic, and there is no chaos analogous to that of the classical 
map: computing the wave function for long times is not more difficult than for short 
times. (There still is hypersensitivity to small perturbations [7], but this is not "chaos" 
in the classical sense.) 

We can also consider quantum systems which are not restricted to a particular value 
of j. Their Hilbert space is spanned by the spherical harmonics Y^ m (^,0), with angles 9 
and related to the cartesian coordinates in Eq. (|l|) in the usual way: x = sin8 cos0, 
etc. The unitary evolution is still generated by Eq. (|3D, where it is now understood that 
the various Jk have to be written as block-diagonal matrices, with blocks of order (2j + 1). 
In that case, it is also possible to write the wave function as 

^(M) = ECWT(M), (4) 

jm 

but its components with different j are never mixed by the U matrix. They evolve 
independently of each other. There still is no chaos: the cost of computing the final 
state does not increase with the number of steps. We thus see that the quantum kicked 
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top is not a faithful replica of the classical T&T map. The reason for this failure of the 
correspondence principle is explained below, and an alternative quantization method will 
be proposed. 

Consider the Hamiltonian which generates the unitary map (^): 

H = aJ 2 /2J + (nJy/2) X) ( 5 ) 

n 

where the unit of time is the duration of one step. In this Hamiltonian, the twist is 
continuous, and the rotation proceeds by kicks (the opposite choice is also possible). If ip 
is represented by a set of Cj m coefficients as in Eq. (|), the Jk in (|) are block-diagonal, 
as explained above. On the other hand we may as well use the "coordinate basis" and 
directly write ip(0, <p) without expanding into spherical harmonics, which is closer to the 
spirit of the classical T&T map. We then have J z = —ihd/d<f), and more complicated 
differential operators for J x and J y [8]. In particular, 



J = —ih 



d ( . „ d\ 1 d 2 11/2 
— sm 9 — + 



sin^ 89 \ 09 sin^fl 



(6) 



When we were in a finite-dimensional Hilbert space, it was natural (indeed unavoidable) 
to write the twist generator as J 2 /2J. However, for arbitray ip(9,<j)), a more natural 
expression for the twist generator is J z cos 9. Let us therefore replace @ by 

H = aJ z cos9 + (7rJ y /2) ]T<5(t-n). (7) 

n 

This is a Hermitian operator which does not commute with J 2 , so that J 2 is not conserved. 
(In other words, J z cos 9 has nonvanishing matrix elements between Yp with different j, 
and cannot be written solely in terms of the Jk matrices.) Instead of (§) we now have 

U = exp(— iirJy/2h) exp(— iaJ z cos 9/h). (8) 

The crucial difference is that the spectra of the new H and U include continuous parts, 
and therefore permit the existence of true chaos. This can be shown as follows. Returning 
to the classical T&T map (1), let us consider, instead of individual points, a mass density 
p(9, 4>) spread on the unit sphere. We may even consider complex densities, if we wish. 
Let us further assume that the infinitesimal mass, p(9, <ft)dQ, attached to the area element 
dQ = sin# d9 d<f>, is conserved by the T&T map. It thus behaves as an incompressible 
fluid, or as a Liouville density would do while moving in phase space. 
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As the map is area preserving, namely d£V = dQ, we have, at each step, 



(9) 



This implies that the map is unitary: for any two densities p\ and p2, 



I 



J 'px{6,(f>y P2 {e,<t>) dn 



(10) 



Such a unitarity property was proved long ago by Koopman [9] for Liouville densities in 
phase space. There is a large body of knowledge on the ergodicity and mixing properties 
of Liouville functions [10]: In the generic non-integrable case, the Liouvillian has a con- 
tinuous spectrum, in which an infinite number of discrete lines may be embedded. If the 
dynamical system has a finite measure, the unitary operator U has at least one eigenvalue 
equal to 1, corresponding to equilibrium. Moreover, it can be proved [10] that, if the 
system is ergodic, but not mixing, that eigenvalue is nondegenerate, and all the other 
eigenvalues of U form a subgroup of the circle group. On the other hand, for a mixing 
system, which also has a single nondegenerate eigenvalue 1, the rest of the spectrum is 
absolutely continuous. A generic dynamical system may have some regions of phase space 
which are subject to mixing, others which are only ergodic, and still others which are not 
even ergodic. Such a system is called "decomposable" [11]. In that case, the spectrum of 
U is continuous, with an infinite number of discrete lines embedded in it. In particular, 
the eigenvalue 1, and possibly others, are degenerate. 

Now, it is easily seen that the quantum wave function ip(6, <fi) behaves, under the 
unitary transformation U in Eq. (|8]), exactly as the fictitious incompressible mass density 
described above. This is obvious for the rotation operator, exp(—iTrJ y /2h), which performs 
a rigid rotation of ip around the |/-axis. For the twist, exp(— iaJ z cos9/ti), we have 



That is, the wave function tp moves on the surface of the unit sphere exactly as the 
classical points did, and it remains normalized, by virtue of the unitarity of the mapping. 
The essential difference between this new quantum system and the former kicked top is 
that the new Hamiltonian (which formally behaves as a classical Liouville operator) has 
a partly continuous spectrum. Therefore the evolution of ij;(9, 0) cannot be represented, 
even approximately, by a finite number of terms. 

Let us examine this evolution more closely. The rotational part of U is represented, in 
the Y™" basis, by the familiar orthogonal matrix U^ m ,, namely a block-diagonal matrix 
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with blocks of size (2j + 1) [12, 13]. In the same basis, the twist (|Tl] ) is also represented 
by a unitary transformation [14], 

Cjm — > C' jm = ^2 Clm, (12) 

I 

where 

Uj^ = J Yf 1 *^, 0) Yr(6, 0) e - macose dVl, (13) 

is a band matrix. This transformation leaves m invariant, but introduces all the j with 
j > \m\ (with exponentially small coefficients for large j). 

It is now possible to give a quantitative measure for the complexity of the wave function. 
Its information entropy [15], 

S = -E \C jm \ 2 \og\C jm \\ (14) 

jm 

has the intuitive meaning that e s roughly indicates the number of basis vectors that are 
appreciably involved in the expansion of ip into spherical harmonics. This "entropy" of 
course depends of the choice of the basis (namely, spherical harmonics). However, it 
can be shown [14] that when S becomes large, its value is asymptotically independent of 
the choice of the basis, provided that the transformation between different bases is given 
by a band matrix (the latter property holds for any two bases whose definition is not 
algorithmically complex) . 

Finally, let us illustrate the above considerations by a numerical example. A plot of S 
versus the number of steps is shown in Fig. 1, for a = 3 and for two initial states, ip±, 
with 

C r 11 = ±C 1) _ 1 = l/A (15) 

respectively, and all other C\ m = 0. It is seen that S increases roughly linearly (with 
a small negative second derivative), so that, as the wave function evolves and becomes 
more and more complicated, the effective number of components needed for representing 
it increases about exponentially with time. (The calculation of the last step in Fig. 1 
involved all the matrices with j, I, m, up to 500. The next step would have exceeded the 
capacity of my workstation, or entailed a severe loss of accuracy.) 

Some of the ideas presented here originated during a visit to the International Center 
of Theoretical Physics, Trieste, whose hospitality is gratefully acknowledged. I had inter- 
esting discussions with Felix Izrailev and Dima Shepelyansky. This work was supported 
by the Gerard Swope Fund, and the Fund for Encouragement of Research. 
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Caption of figure 

FIG. 1. Growth of S for the two initial states in Eq. (|15|). 
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